Vortex wandering in a forest of splayed columnar defects 
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We investigate the scaling properties of single flux lines in a random pinning landscape consisting of 
splayed columnar defects. Such correlated defects can be injected into Type II superconductors by 
inducing nuclear fission or via direct heavy ion irradiation. The result is often very efficient pinning 
of the vortices which gives, e.g., a strongly enhanced critical current. The wandering exponent £ 
and the free energy exponent ui of a single flux line in such a disordered environment are obtained 
analytically from scaling arguments combined with extreme-value statistics. In contrast to the 
case of point disorder, where these exponents are universal, we find a dependence of the exponents 
on details in the probability distribution of the low lying energies of the columnar defects. The 
analytical results show excellent agreement with numerical transfer matrix calculations in two and 
three dimensions. 

PACS numbers: 74.60.Ge (Flux pinning), 02.50.-r (Prob. theory), 74.60.-w (Type-II sup.) 



I. INTRODUCTION 

Pinning of vortex lines to defects in superconductors 
plays an eprtremely important role in determining their 
properties u Motion of flux lines in response to an applied 
current induces a voltage and hence dissipation. Only 
if the flux lines remain pinned will the linear resistance 
be truly zero and the material exhibit superconductivity. 
By introducing controlled disorder it is possible to signif- 
icantly improve the critical currents, temperatures, and 
fields. A particularly efficient strategy is to bombard the 
sample with heavy ions to produce linear damage tracks 
which form optimal pinning centers for the vortex lines .□ 
When these columnar defects are parallel the vortices 
will, below a critical temperature, become localized in 
a Bose glass phase that replaces the Abrikosov vortex 
lattice of the clean systems Further improvements were 
suggested to occur for splayed columnar defects,cl where 
the angle mismatch between neighboring columns should 
lead to increased barriers for variable range hopping and 
therefore an even stronger enhancement of the critical 
current in the splay glass phase. Even more important is 
perhaps the suppression of vortex motion due to forced 
entanglement that is created by the splay. 

There are several ways in which splayed configurations 
of columnar defects can be produced. A narrow approx- 
imately Gaussian angle distribution can be produced by 
placing a thin metal foil in front of the sample during 
irradiation to defocus the beam of heavy ions as they en- 
ter the sample. Another interesting approach is to induce 
fission in some of the nuclei in the material. The fission 
products will then go apart in opposite directions cre- 
ating damage tracks with a more isotropic distribution 
of angles. Lastly, heavy ion irradiation can be applied 
at different discrete angles to create several families of 
parallel columnar defects. In this way it has been possi- 
ble to increase critical currents by more than an order of 
magnitude in many cases. □ 

The theoretical understanding of the Bose glass, occur- 



ring for parallel columnar defects, rests to a large extent 
on the mapping of the vortex line problem to a quantum 
mechanical system of (2+l)-dimcnsional bosons in a ran- 
dom potential, with imaginary time playing the role of 
the z-axis parallel to the external magnetic field. With 
randomly splayed defects, however, the disorder devel- 
ops long range correlations in space and "time" , and this 
analogy is less useful. In fact, splay disorder leads to 
logarithmically divergent phase fluctuations in the boson 
order parameters The properties of individual flux lines 
in such disordered environments are just beginning to be 
investigated.!!! 

In this paper we focus on the scaling properties of a 
single flux line in a sample with many randomly tilted 
columnar defects. We consider the flux line as it enters 
the superconducting sample on one side at some arbitrary 
but fixed position (taken to be the origin), and leaves 
on the opposite. Physically, this fixing of the starting 
position could be due to, e.g., surface pinning as in the 
experimental setup of Fig. |l| below. The wandering and 
energetics of the flux line can be characterized by two 
exponents defined by the following relations: 



(r 2 (*)) - 
AF 2 (z) 



,2C 



(la) 
(lb) 



where r is the transverse position, z is the distance tra- 
versed parallel to the magnetic field, F is the free energy, 
AF = F — F, the overbar denotes a disorder average, 
and (■ • •) a thermal average. 

These problems, usually referred to as directed poly- 
mers in random media ODPRM) , have received much at- 
tention in recent years ,u with interesting connections to 
many other very different physical systems. The prob- 
lem can, e.g., be mapped to the noisy Burgers equa- 
tion describing a randomly stirred turbulent fluid, or 
equivalently, the Kardar-Parisi-Zhang (KPZ) model of 
surface growth. Without disorder, a flux line subject 
only to thermal fluctuations will make a diffusion-like 
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random walk as it wanders from the bottom to the top 
of the sample with C — 5 an d u) = 0. Point disorder 
tends to increase the transverse wandering of the flux 
line leading to a non-trivial universal exponent £ > i. 
(In (1+1) dimensions the exponents are known to be ex- 
actly C = 2/3. and lu — 1/3, in (2+1) we have £ ps 5/8 
and u> w 1/4.0) In this case statistical tilt symmetry en- 
sures that ui — 2£ — 1 for arbitrary d. The fact that 
u> > suggests that energy barriers become arbitrarily 
large for large z, so that the system is governed by a 
zero temperature fixed point, where thermal fluctuations 
are irrelevant on long length scales. Parallel columnar 
disorder, on the other hand, tends to localize the flux 
line, i.e., to reduce its transverse fluctuations. However, 
in search of ever lower pinnins-,energies it will still wan- 
der with non-trivial exponents.Q At zero temperature the 
exponents depend on the low energy details in the dis- 
order distribution, e.g., £ = lo = v j{d + v) for bounded 
distributions of the form Eq. (ji|) below, while at finite 
temperature the wandering is slightly sub-ballistic, with 
|r(z)| ~ z/(lnz) 7 , 7 = 1 + 2/d, .where d is the number 
of dimensions perpendicular to zn For splayed columnar 
defects one might expect an enhancement of this wander- 
ing since the flux lines try to follow the randomly tilted 
defects. 

The path a flux lines takes inside a superconductor is 
not easy to measure experimentally. What can be mea- 
sured is the positions at which the flux line enter and 
leave the superconductor, e^g., using double-sided deco- 
ration experiments. In Fig. |f]we propose an experimental 
setup which would allow the endpoints r(z) of the flux 
line to be measured for several different disorder realiza- 
tions of varying height z, using a single piece of super- 
conductor. Short parallel columnar defects on one side 
of the sample provide a set of known starting points for 
the vortices as they enter the forest of splayed columnar 
defects and traverse the sample. At low magnetic field 
it should be possible to approximately matciL-entry and 
exit points via double-sided flux decorations .113 The vari- 
able length parallel defects simulate the effect of varying 
sample thickness. 

Apart from the intrinsic interest in the behavior of 
a single flux line in a disordered superconductor, the 
physics of a single vortex in the presence of splay is 
relevant for understanding properties at low magnetic 
fields, where the flux density is low. It affects, e.g., 
the constitutive relation B(H) just above the -ff c i-line, 
where a balance between magnetic field, pinning energy, 
and repulsive interactions gives B ~ (H — i? c i) > with 
(3 = dQ/ (1— w).lllr0 This relation was recently measured in 
experiments in a (1+r-U dimensional geometry for samples 
with point disorder ,113 confirming the prediction (3=1. 
Furthermore, the non-linear current-voltage (IV) charac- 
teristics at low temperatures and currents is determined 
by the depinning of single flux lines from the defects. If 
we assume that the free energy barriers scale in the same 
way as the sample-to-sample free energy fluctuations, an 
excitation of length z and width £ ~ zS would cost pin- 




FIG. 1. Experimental setup that would allow measure- 
ment of the wandering exponent Parallel columnar defects 
of varying depth are introduced from the bottom of the sam- 
ple in order to provide a strong pinning mechanism to fix the 
entry positions of the flux lines. Splayed columnar defects 
are introduced from the other side, and the relative trans- 
verse distance between the starting point and the exit point 
of the flux line could then be measured, e.g., in double sided 
decoration experiments, or using magneto-optical imaging. A 
similar setup could also be used to investigate the vortex wan- 
dering in point disordered samples. 

ning energy AE ~ z u but gain J£z due to the Lorenz 
force from the current density J, so that the free energy 
barrier to overcome is AF(J, z) = k\z^ — kiJz 1+t * . Op- 
timizing with respect to z gives AF(J) = AJ~^, where 
/i = w/(e + l — a;). Thermally activated creep over the 
free energy barriers then gives rise to an electric field 
£ ~ exp[-(J (T)/Jn with J (T) = (A/T) 1 ^. Fi- 
nally, the study of single flux lines also allows estimates of 
entanglement lengths of interacting vortices, which may 
shed light on the still mysterious properties of the splay 
glass phase. 

While the single vortex pinning from point disorder 
and parallel columnar defects has received much atten- 
tion, the case with splayed columnar disorder has only 
been addressed quite recentlyj] and only for the spe- 
cial case of a "nearly" isotropic distribution of splay 
angles. There it was argued that, for very wide angle 
distributions, the problem would reduce to a flux line 
in a disorder landscape with long range spatial trans- 
verse correlations but no correlations along the mag- 
netic field. Explicitly, the Fourier transform of the 
disorder correlator for general splayed configurations 
is A(k±,k z ) = A f(k z /v k±)/v k±, for small k ±1 k Zl 
where vq is the characteristic width of the distribution 
of slopes. "Nearly isotropic" splay refers to the limit 
vq — ► 00 (keeping Aq/vq fixed), giving A(k±,k z ) oc 
1/fex, while a truly isotropic distribution would give 
A(k±,k z ) cx 1/ \/k'j_ + fef . Neglect of the fc z -dependence 
is at least plausible in the case £ < f, since the width 
grows like v ~ [length] under renormalization and 
therefore would diverge to infinity on large length scales. 
Subject to these assumptions, the flux line problem can 
then be mapped to ajaoisy Burgers equation with sp.a=. 
tially correlated noise,Q for which the universal resultsE3 
C = 3/(3 + d) and ui = (3 — d)/(3 + d) have been con- 
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jectured. In this approach, however, higher order corre- 
lation functions of the disorder are neglected. 

Here we consider splayed columnar defects for general 
angle distributions using a real space approach, similar in 
spirit to Ref. || for parallel columnar defects. This leads 
to a rather different physical picture, where the impor- 
tant feature is that the flux line keeps following colum- 
nar pins with lower and lower energy (and smaller and 
smaller slopes). We find new exponents that depend on 
the tail of the distribution of the energies of the columns. 
The sensitivity to low energies suggests that higher cor- 
relation functions of the disorder, which were missing in 
the previous treatment, may be important. The predic- 
tions for C and uj of the different approaches are compared 
in Table Q. We have been able to confirm the results of 
the real space approach using numerical transfer matrix 
calculations in (1+1) and (2+1) dimensions. 

In Sec. we use a combination of scaling arguments 
and extremal statistics to derive analytic expressions for 
the wandering and energy exponents. We also discuss 
crossover effects and other complications, such as inho- 
mogeneities of the columnar defects (fragmentation) , and 
the influence of temperature. Section fin describes our 
numerical work. We introduce a model well suited for 
transfer matrix calculations of superconductors contain- 
ing splayed columnar defects, discuss finite size effects, 
and present our results. In Sec. [jy] we review some open 
questions. 



II. ANALYTIC ARGUMENTS 

A. The model 

The Hamiltonian of the flux line consists of an elastic 
energy and an interaction with the disorder, 



H 




+ U(r(z),z) 



(2) 



where r is the transverse d-dimensional position, z 
the distance parallel to the magnetic field, and 7 is 
the line tension. Splayed columnar defects are dis- 
tributed throughout the system, giving a disorder poten- 
tial U(r,z) — Ui(r — Viz — r?), with random slopes 
Vj and random potential wells U%. The r°'s are the posi- 
tions where the columns cross the plane z — 0. We are 
interested in the low-temperature regime where the flux 
line is tightly bound to the defects. The binding eneig 
per length of a flux line pinned to a column is thercr 
£i = fvf + m, where u t « U t (0) + cT 2 / (2-fb 2 ) , b t being 
the radius of the defects, and c a constant of order unity. 

The probability distribution P(e) of the binding en- 
ergies depends on details of the injection process in the 
material. Below we will need information about the low 
energy tails of this distribution. Given the probability 
distributions P u (u) of the depths and P«(v) of the slopes, 



which we assume to be independent, we get for the prob- 
ability of having total energy per length e the expression 



P(e) = I du I d d vP u {u)P v {v)6{e-u-]-v 2 ) 



(3) 



In any realistic case the potential distribution will be 
bounded, and for convenience we will use 



— , if < u < A 



PJu) = { A" 

1 0, otherwise. 



(4) 



A uniform bounded distribution is perhaps the most nat- 
ural choice and corresponds to v — 1, while the limit 
v — > gives a delta function, i.e., all columns have the 
same depth. In the general case v > 0. We have further 
chosen the energy scale such that zero corresponds to the 
deepest possible pin. We assume a continuous symmetric 
distribution, P„(v), of the slope of the columnar defects, 
with a finite value for v — > and a characteristic width 
i>o . Low energies turn out to be the most important ones, 
and we find 



< e < Zv? 



! 



A, 



i.e., a power law behavior P(e) ~ £ Q 1 with 

d 



(5) 



(6) 



Under the assumptions made above this relation should 
hold in the tail even in the presence of weak correlations 
between P u {u) and P v (v). 

A complication may arise from the fact that the defects 
are not necessarily uniform along their extension. Since 
the heavy ions creating the damage tracks lose energy 
to the material as they travel through it, there may be 
a systematic dependence of the diameter and angle on 
z. The tracks formed by low-energy ions are furthermore 
often fragmented. We will neglect these effects for now 
and return to them in the discussion below. With a good 
model of these processes it should be possible to extend 
the arguments given below to take into account these 
facts. In practice P(e) could have several regimes with 
an approximate power law behavior, leading to crossovers 
between different scaling regimes before the asymptotic 
result determined by Eq. (pf), is reached (see below). 



B. Scaling theory from extremal statistics 

As a flux line travels through the sample at low tem- 
perature it seeks out pins with lower and lower energy. At 
a given distance z the flux line has had a chance to find 
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the optimal columnar pin within a region of size l d x z, 
with £ ~ z*>. In this region there are on average N = pi 
pins, where p is the areal density perpendicular to z. For 
a given length scale i the physics will be dominated by 
the pin with lowest energy among these N pins. The 
typical energy per unit length of this optimal col- 
umn can_then be estimated using arguments of extremal 
statistical if N is large enough, i.e., for large z. Upon 
introducing the probability that a pin has an energy less 
than e, 



F(e) 



P(e')de' oc 



A" ( 7Wo 2 /2) 



d/2 • 



(7) 



we have 

Prob(e opt > e) = (1 - F(e)) N 



-NF(e 



for N large. 

(8) 



For a typical energy of the optimal column we have 
NF{<%>) ~ 1, i.e., ' 



typ 
fc opt 



N~ 



(9) 



The total energy of the whole optimal path will then scale 
as 



typ 



fc opt 



(z')dz' 



Thus we have obtained the relation 



uj = 1 — (d/a 



(10) 



(11) 



between the energy exponent and the wandering expo- 
nent. Note the dependence (through a and Eq. (0)) on 
the detailed behavior of P u {u) at low energies. 

Now we will turn to the determination of £ itself. As- 
suming that the flux line is mostly localized along the 
columns with lowest energy, we only need to know how 
the slope v op t of the optimal pin encountered decreases 
with z. The probability of v opt is, for large N, given by 

Pv opt (v op t)d d v opt = 

N [ deP u (e - lv 2 opt )P v (v opt )e- NF ^d d v opt . (12) 



A simple calculation shows that ~ J 2e^ p \/^ ~ 

N -i/2 a _ z -(d/2 a for £ typ < 7W 2/2, and for the typi- 
cal value of the endpoint r(z) of the flux line at distance 
z we then obtain 



(13) 



ThusC = l/(l + d/2a) = 2a/(2a + d), and via Eq. (jllj), 
= (2a — d)/(2a + d). From Eq. (^) we then obtain our 
final result, 



TABLE I. Summary of wandering and energy exponents 
for some selected values of v and for other kinds of disorder. 
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3/4 
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(14a) 
(14b) 



Remarkably, these exponents satisfy for all d and v the 
relation u> = 2£ — 1 valid for systems obeying statistical 
tilt symmetryja even though strictly speaking, there is no 
exact tilt symmetry in the present case. 

Since < v < oo, we always have \ < ( < 1 and 
< u> < 1. For a uniform bounded disorder distribu- 
tion (y = 1) we get C = 3/4, uj = 1/2 in (1+1)D and 
C = 2/3, uj = 1/3 in (2+l)D, leading to an enhanced 
wandering compared to point disorder. If all splayed col- 
umn well depths are identical, v — > giving £ = | and 
uj = 0. This means that the T = fixed point becomes 
marginal and thermal fluctuations (as well as point disor- 
der) could become important. For d — > oo the exponents 
smoothly attain their thermal values £ = ^, uj = 0. If 
P u (u) is replaced by a Gaussian distribution (thus allow- 
ing arbitrarily low pinning energies) we find £ = uj = 1, 
described by the v — > oo limit, but in this case there 
may also be logarithmic corrections. A distribution with 
power-law tails, P u (u) ~ u < — 1, would give 

( = 1 and ui = 1 + d/v. In Table || we summarize the 
exponents found for splayed columnar defects and com- 
pare them with the exponents for other kinds of disorder. 
Comparing with the results found in Ref. [?] for a disor- 
der distribution with long range transverse correlations, 
A(fc) ~ l/fc±, the exponents agree for the special case 
d = 1, v = 1, but disagree in general. 

Another interesting quantity to look for is how often 
the flux line changes direction. Let us define the total 
number of intercolumn jumps the flux line has made up 
till z as n z . The turning probability in an infinitesimal 
interval [z, z + dz] is given by the number of new columns 
explored times the probability of any of them being lower 
in energy than the current energy, dn z ~ F{e^ pt )dN ~ 



dN/N, which gives 



\nN(z) - lnz. 



(15) 



From this follows that the spacing along z between two 
turns increases like Az ~ z on average. 
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C. Crossover effects and scaling regimes 

The results above hold asymptotically for large z and 
depend only on the parameter v in the pinning energy 
distribution. Realistic pinning energy distributions can 
be expected to be more complicated than the simple form 
Eq. (y) used above and this may lead to various crossover 
effects between different scaling regimes, depending on 
the energy scale probed at a given length scale. This 
may effectively lead to a z-dependent parameter v, and 

and hence the 



may also disrupt the relation 
relation lu = 2Q — 1 over certain length scales, although 
asymptotically the behavior should be that of the previ- 
ous section. 

For the bounded distributions used above, Eq. (|[), the 
asymptotic behavior sets in when < min {A, 7Dq/2}, 
or, using 



2„2\ d / 2 



(16) 



where N(z) is the number of columns explored up till z. 
Finally, we expect behavior typical of parallel columnar 



defects when z <C aj_/v , where a± = p 
this crossover scale diverges as expected. 



-l/d 



As v -> 



D. Fragmented columnar defects 

In some cases the columnar defects created by the ion 
irradiation (or fission fragments) are highly nonuniform 
along their axis (See Fig. ||).E£| This fragmentation de- 
pends on, e.g., the kinetic energy and the type of ions 
used. In this subsection we discuss the influence of this 
fragmentation on the flux line wandering. 

W e firs t consider the stability of the results derived in 
Sec. [IB. With fragmentation the energy of a flux line 



pinned to a columnar defect will be subject to two kinds 
of disorder. As before there is a random z-independent 
part constant along the columns, given by Eq. (J4j) . We 
assume that fragmentation superimposes a random short 
range correlated energy modulation along the columns 
with zero mean. 

Let us now consider a segment of a flux line pinned 
to a defect over a length Az 3> a, where a is the mi- 
croscopic size of the fragments. The fluctuation in the 
average energy per unit length then gets a contribution 
~ ef = f^/a/Az from the fragmentation, where / is 
the root mean square strength of the disorder. Frag- 
mentation should be relevant only if this energy exceeds 
the fluctuations in the constant part of the energy. The 
density of columns with energy less than this energy is 
p < = pF(ef), and their mean transverse distance t m 
{p<)~ l ' d . The mean longitudinal distance between cross- 
ings (or close encounters) of these low energy columns is 
then Az* ss l/v*, where the v* is the maximum typical 





FIG. 2. Illustration of a couple of fragmented columnar 
defects (left) and a flux line which fluctuates back and forth 
between two crossing columnar defects at finite temperature 
(right). The entropic contribution to the free energy from the 
crossing points, which appear randomly along the defects, act 
very much like fragmentation. 



slope, v* ~ y / 2ef/"f ~ y/2f/~f (a/Az) 1 / 4 . Introducing 
the mean transverse spacing of defects ai = p~ x l d we 
find 



A: 



a± ThL^/d f Az 



vo f 



1 + v/d 



(17) 



Fragmentation can only be important if Az* < Az, i.e., 
on large length scales if v < d, and otherwise on short 
length scales, with the crossover scale given by the re- 
quirement Az* « Az, 



Az, 



crossover 



voa f 



l—v/d 



l+v/d 



(18) 



What happens when the fragmentation dominates, i.e., 
on long length scales for v < d or on short length scales 
for v > dl Although the disorder along the columns 
is correlated only over the microscopic scale a of the 
fragments the problem still has long ranged correlations 
in that the disorder is restricted to the tracks of the 
columns. There are (at least) two possible scenarios. 
If the point-like contribution to the disorder from many 
columnar fragments is strong enough to completely de- 
localize the flux line, this would presumably lead to a 
crossover to the physics dominated by point disorder, 
with a corresponding change in the critical exponents. 
However, if the fragmentation is weak, the flux line could 
remain pinned to the best fragmented column available 
in the region explored, occasionally making excursions 
to take advantage of favorable fluctuations in the dis- 
order. In Ref. 16 the depinning of a flux line from a 
single columnar defect in the presence of point disorder 
was studied. These authors found a transition from a 
localized to a delocalized state for d > 1 as a function 
of disorder strength, while the case d = 1 was marginal 
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with the flux line pinned by an arbitrarily weak attractive 
column. In analogy with these results we do not expect 
the first scenario to be realized for d = 1, whereas in 
higher dimensions cither scenario is a priori possible and 
a phase transition between the two cannot be excluded. 

We now take a closer look at the second scenario of 
weak fragmentation. The flux line will then effectively see 
splayed columnar defects with a Gaussian distribution of 
pinning energies with variance ~ f 2 a/ Az, that explicitly 
depends on the distance z (recall that Az ~ zL Repeat- 
ing the extremal statistics calculation of Sec. [II B| for t his 
choice of pinning energies leads to ~ f^/a/Az up 
to logarithmic corrections. The typical optimal slope is 

«opt ~ ^2e*pt/7 again ignoring logarithms. The total 
energy of the flux line then scales as E typ (z) ~ z 1 / 2 , and 
the endpoint position as r typ (z) ~ z 3 / 4 , so that 

C = 3/4 (19a) 
a; = 1/2, (19b) 

when weak fragmentation dominates. 

E. Finite temperature 

We have assumed thus far that temperature was low 
enough to be ignored and the problem could be ana- 
lyzed solely in terms of energetics. As the tempera- 
ture is increased entropy becomes more important. At 
relatively low temperatures the main contribution will 
come from regions where two or more columns come 
closer than It » T/y / 2~fAU, where AU is the energy 
difference between a columnar defect and the surround- 
ing material. The flux line will then gain entropy by 
fluctuating back and forth between these columnso (See 
Fig. ^]) . Since these crossing points occur randomly along 
the splayed columns the situation becomes very similar 
to that of fragmented columnar defects at zero temper- 
ature. Therefore, we expect the argument of the previ- 
ous subsection to apply, i.e., temperature will be impor- 
tant whenever fragmentation would be important, and 
the wandering and energy exponents should be the same 
as for that case. If the temperature is increased even fur- 
ther the flux line will no longer be pinned to individual 
defects, but rather collections of them. In this case the 
situation becomes far more complicated, and beyond the 
scope of the present paper. 

III. TRANSFER MATRIX CALCULATIONS 

To check the analytic arguments in Sec. ||, we have 
performed numerical transfer matrix calculations. Usu- 
ally transfer matrix calculations of directed polymers are 
performed on a model defined on a regular lattice, but 
for the present problem it is difficult to account for the 
varying splay angles using this approach. Instead we have 



developed a model on an irregular random lattice defined 
by a particular configuration of randomly splayed colum- 
nar defects. 



A. Model 

The disorder landscape consists of N splayed columnar 
defects embedded randomly in a system of size L d x L z . 
The transverse positions of the columnar defects are 
given by 

Ti(z) = r s + v 2 z, (20) 

where the tj are distributed uniformly in the plane z = 0, 
and the Vj's are chosen randomly with some distribu- 
tion P(v)d v. We also assign a random energy cost 
Ui 6 [P>A1 to each column, with probability P u (u)du 
from Eq. (|4j). Periodic boundary conditions are employed 
so that we always have Ti(z) € [-L/2, +L/2] d . In addi- 
tion we use a discretization in the z-direction with lattice 
constant a z = 1. 

A flux line enters the system at Ro = (0, 0) and leaves 
at some R = (r, L z ). The model is defined by restricting 
the positions of the flux line to the columnar defects. The 
space between pins is excluded, except to allow jumps 
from one column to another (see Fig. ||). The flux line 
can then be parameterized by the sequence {iz} z =o °f 
columns visited as it traverses the sample from bottom to 
top. The actual path a vortex takes is given by r^ z (z) = 
rf + Vj z, and the total energy of the flux line is given 
by 

L z -1 

H = AEi M+1 , it (z), (21) 

2=0 

where AEy (z) is the energy of the (straight) line segment 
between r.;(z + 1) and Yj(z). We take 

AE^z) = 1 [r 4 (z + 1) - v 3 {z)f + ui, (22) 

where the first term is the elastic energy (the distance ap- 
pearing in brackets is understood to be the shortest dis- 
tance using the periodic boundary conditions) and the 
second is a random site energy of column i, consistent 
with the continuum model Eq. (Q). More elaborate forms 
could be used but will presumably not change any uni- 
versal properties. Fragmentation of the columnar defects 
is modeled below by adding an uncorrelated uniformly 
distributed term fi(z) € [0,/] to Eq. (||). 

B. Recursion relations 

The partition function for the flux line now obeys the 
recursion relation 



G 




FIG. 3. Discretization used for transfer matrix calcula- 
tions in (1+1) dimensions. At each step the flux line has 
the possibility to jump to anyone of the other columnar pins 
at the next constant-,? section. A similar picture obtains in 
(2+1) dimensions. 

JV 

Zi(z + 1) =J2 e ~ 0AEij(z)z i( z )> (23) 

3=1 

where (3 = 1/T is the inverse temperature, and we use 
units such that fee = 1 ■ At zero temperature this reduces 
to an optimization problem for the total energy of the 
path: 

Ei(z + 1) = min{^-(z) + AE^z)} (24) 
j 

The sum in Eq. ( p3[ ) and minimization in Eq. ( p4| ) are 
over all possible columns, thus jumps between columns 
are not limited to just nearest neighbors, but instead re- 
stricted by the elastic line tension. These relations are 
easily iterated on a computer, and averages are then cal- 
culated from 




FIG. 4. Example of the optimal paths with variable end- 
points in a particular disorder realization in (1+1) dimensions. 
The light gray lines represent the randomly splayed columns 
with random energies given by Eq. (^) with v = 1. The solid 
lines are the optimal paths if the endpoint is restricted to a 
given position. The thickest line is the optimal path among 
all of them, i.e., without any restrictions. The upper panel 
shows the corresponding energies as a function of endpoint 
position. 



(0«> 



(25) 



The free energy is given by F(z) = —T\a.Z{z), where 
Z(z) = '^2 i Zi(z). (At zero temperature the free energy 
is of course equal to the energy.) The overbar denotes 
the quenched average over different disorder realizations. 



hand sides become independent of L. In the opposite 
limit of large arguments C(x) — > d/Y2, and f(x) ~ a; 2 .Ell 
The crossover happens at some value of x = x r 



i.e. 



for z = x c L 1 ^. Thus, when (r 2 )/L 2 and AF 2 /L 2x are 
plotted against z/L 1 ^, the data for different L should 
collapse onto a single curve (at least for large enough z). 



C. Finite size scaling 

In a system of finite transverse size L, we expect that 
the mean square fluctuations of the endpoint will obey 

Fi^ 2 CWi 1/( ), (26) 

while the free energy fluctuations satisfy 

~KFHz) = L^f{z/L 1 ^), (27) 

where \ — W /C- For small arguments the functions C 
and / should become power-laws, such that the right 



D. Results and discussion 

The recursion relations Eq. (^3|) and ( pi| ) were solved 
numerically and averaged over 2000-10000 disorder real- 
izations to get small error bars. Figure ^| shows an exam- 
ple of the optimal path (corresponding to T = 0) of a flux 
line for a particular disorder realization. The number of 
defects were taken to be TV — L d so that the mean spacing 
between defects is unity. The line tension is set to unity, 
equivalent to measuring all energies in units of 7. We 
further choose the distribution of tilts P(v) uniform in 
[— vq, +vo\ d . Apart from the shape of the tail of the pin- 
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ning potential, given by the exponent v, the asymptotic 
behavior on long length scales should be independent of 
A and vq. However, these values will influence the length 
scales where the system crosses over to the asymptotic 
regime. Since the system sizes possible to study numer- 
ically are limited, care must be taken in choosing the 
values of these parameters to reduce crossover effects, 
otherwise inaccurate results can easily be obtained. The 
analysis of Sec. [IC suggests that the crossover length 
scales are minimized when « A. With this in 

mind we usually put vq = 1/2, A = 0.125 and sometimes 
w = 1/4, A = 0.03. 

The inset of Fig. 0(a) shows the mean square positional 
fluctuations as a function of z for d = 1 and v — 1 at zero 
temperature. Note that the curves saturate for large z, 
consistent with expectations from Eq. (^6|). The main 
part of the figure shows a scaling collapse of the same 
data using the theoretical value of £ = 3/4. In Fig. ||(b) 
we show a scaling collapse of the energy fluctuations using 
ijj = 1/2. Here there are clearly visible corrections to 
scaling for small z, and no power-law behavior is observed 
for the smallest system sizes. However, as L is increased, 
the range over which scaling occurs increases. In Fig. 0(a) 
and 0(b) we show similar results for d = 2 and v = 1, 
where Eq. ( fu] ) predicts C = 2/3 and ui — 1/3, and again 
the agreement is excellent. 

To test the dependence of the exponents on the tail 
of the pinning energy distribution we calculate the mean 
square positional fluctuations [Fig. 0(a)] and the energy 
fluctuations [Fig. 0(b)] for several different values of v. 
For large z the data is accurately described by power 
laws with exponents given by Eq. ([14]) . It is also possible 
to collapse the data for different system sizes in each of 
these cases in a similar way to Fig. g. 

In Fig. 0we show results for fragmented splayed colum- 
nar defects. In (a) we plot the positional fluctuations for 
d = 1, v = 0, for increasing strengths of the fragmenta- 
tion /. Since v < d fragmentation should be relevant, 
and already for small values of / we see deviations from 
the f — result, where £ = 1/2. As f is increased further 
£ increases and stabilizes to a larger value. For large val- 
ues off a power-law fit gives C ~ 3/4 in agreement with 
Eq. (19). In (b) we plot the same quantity for d = 2. For 
these parameters the scaling crosses over to that of point 
disorder with £ ss 5/8. Thus, fragmentation seems to 
lead to point-like wandering for d — 2, while the results 
for d = 1 are consist ent w ith the second weak fragmenta- 
tion scenario of Sec. II D . For v > d weak fragmentation 
should be irrelevant, and we do find that the exponents 
for homogenous columns, Eq. (|l4|), are much more sta- 
ble for small / than they are for v < d. For large / the 
behavior is still dominated by fragmentation on short 
length scales. An accurate verification of the crossover 
scale, Eq. (|l8|), is very difficult with the limited system 
sizes we have studied, and has not been attempted. 

Finally, we explore the influence of finite temperature 
in Fig. 0. As proposed in Sec. HE the scaling shows 
strong resemblance to that of fragmented columnar de- 



fects, compare Fig. 0. In (1+1) dimensions [Fig. 0(a)] a 
powerlaw fit gives £ « 0.75 — 0.8, which is consistent with 
3/4, the value for fragmented columnar defects. In (2+1) 
dimensions [(b)] we get instead £ « 0.6 — 0.65, consistent 
with the value for point disorder, £ w 5/8. 



IV. FUTURE DIRECTIONS 

Our results indicate that the values of the wandering 
and energy exponents depend sensitively on the low en- 
ergy properties of the pinning energy distribution. In ex- 
periments one might therefore expect these exponents to 
depend on the method of sample preparation. A detailed 
comparison with experiments would require information 
about the pinning energy distribution of the columnar 
defects for the particular samples studied. In case the 
columnar defects are fragmented (and v < d) this de- 
pendence should disappear and universality be restored. 
The prospect of having a phase t ransi tion between the 
two fragmentation scenarios of Sec [I D is interesting and 
deserves further study. 

In this paper we have focused our attention on the 
properties of single flux lines in superconductors with 
splayed columnar defects. At high vortex densities, i.e., 
at high magnetic fields, interactions between vortices can 
no longer be neglected and the problem becomes signif- 
icantly more complicated. It would be interesting, al- 
though not easy, to extend both the analytic arguments 
and the transfer matrix calculations to the case of many 
interacting vortex lines. We do, however, expect that the 
sensitivity to the low energy tails of the pinning energy 
distribution should go away as the vortex density is in- 
creased, since the low energy pins will all be occupied in 
this case. 
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FIG. 5. (a) Scaling collapse according to Eq. (|26|) of mean square fluctuations of the endpoint position of a flux line in (1+1) 
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is observed. The inset shows the unrescaled data, i.e., r 2 vs z. (b) Scaling collapse of the energy fluctuations of the optimal 
path. The dashed line is a power law with exponent 2u> = 1. The inset shows the unrescaled data, i.e., AF 2 vs z. 
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FIG. 6. (a) Scaling collapse of mean square fluctuations of the endpoint position of a flux line in (2+1) dimensions for 
different system sizes at T = 0. A power law behavior with exponent 2£ = 4/3 (dashed line) is observed. The inset shows the 
unrescaled data, i.e., r 2 vs z. (b) Scaling collapse of the energy fluctuations of the optimal path. The dashed line is a power 
law with exponent 2uj = 2/3. The inset shows the unrescaled data, i.e., AF 2 vs z. 
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FIG. 7. (a) Mean square positional fluctuations for several different pinning energy distributions, Eq. (gj), parameterized 
by v. Power-law fits (not shown) to the straight parts of the curves give exponents that agree with Eq. (114) to within 1%. (b) 
Energy fluctuations for different v. 
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FIG. 8. Positional fluctuations for fragmented splayed columnar defects in (a) (1+1) dimensions and (b) (2+1). As the 
strenght of the fragmentation increases the wandering exponent changes from ( = 1/2 to a larger value. The dashed lines are 
power laws with exponents 5 = 3/4 and 5 = 5/8, respectively. 
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we would instead have found {r 2 (z)) 



v 2 z 2 

y opt^ • 



In fact, 



combining this result with the free energy, which goes like 
AF 2 (z) — > Ae 2 pt z 2 , and the finite size scaling relations 
Eqs. ( pqp7j ), gives us an alternative way of arriving at the 
results of Sec. II B, since the L dependence of Ae 2 pt and 



"opt 



can be calculated from extremal statistics. 
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